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Abstract 

We present a calculation of the complete electroweak 0(a) corrections to 
pp — > — > i^uX (£ = e, fi) in the Standard Model of electroweak in- 
teractions, focusing on those corrections which do not contribute in the pole 
approximation. We study in detail the effect of these corrections on the trans- 
verse mass distribution, the TV-width measurement, and the transverse mass 
ratio and cross section ratio of W and Z bosons. 
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I. INTRODUCTION 



The Standard Model of electroweak interactions (SM) so far withstood all experimental 
challenges and is tested as a quantum field theory at the 0.1% level [1]. However, the 
mechanism of mass generation in the SM predicts the existence of a Higgs boson which, so 
far, has eluded direct observation. Direct searches at LEP2 give a 95% confidence-level lower 
bound on the mass of the SM Higgs boson of Mh > 114.4 GeV [2]. Indirect information 
on the mass of the Higgs boson can be extracted from the Mh dependence of radiative 
corrections to the W boson mass. With the present knowledge of the W boson and top 
quark masses [1,3], and the electromagnetic coupling constant, a(M§) [4], the SM Higgs 
boson mass can be indirectly constrained to Mh = 113^2 GeV [1,5] by a global fit to all 
electroweak precision data. Future more precise measurements of the W boson and top 
quark masses are expected to considerably improve the present indirect bound on M H - with 
a precision of 27 MeV for the W boson mass, M W) and 2.7 GeV for the top quark mass, 
which are target values for the expected integrated luminosity of 2 ftr 1 of Run II of the 
Tevatron, Mh can be predicted with an uncertainty of about 35% [6]. In addition, the 
confrontation of a precisely measured W boson mass with the indirect SM prediction from 
a global fit to all electroweak precision data, M w = 80.386 ± 0.023 GeV [1] will provide a 
stringent test of the SM. A precise measurement of the W width, T W) and comparison with 
the SM prediction, will help to further scrutinize the SM. The W width is expected to be 
measured in Run II with a precision of about 25 — 30 MeV from the high transverse mass 
tail when data from both lepton channels and both experiments are combined [7]. It can 
also be determined indirectly from the cross section ratio [8,9], 

Rw/Z ~ a( P p ^z^ e+e-x) ' (1) 

together with the theoretical prediction for the ratio of the total W and Z production 
cross sections, the LEP measurement of the branching ratio B(Z — > [1], and the SM 

prediction for the W — > iv decay width. 

In hadronic collisions, the W mass is usually determined from the transverse mass distri- 
bution of the lepton and neutrino which originate from the W decay, W — > Iv. The experi- 
mental uncertainty on Mw strongly depends on an accurate measurement of the transverse 
momentum of the neutrino which requires that the transverse momentum distribution of the 
W is well understood. In lowest order, the W is produced without any transverse momen- 
tum. Only when QCD corrections are taken into account does the W acquire a non-zero 
transverse momentum, p^ . For a detailed understanding of the p^ distribution, it is nec- 
essary to resum the soft gluon emission terms [10], and to model non-perturbative QCD 
corrections [11]. In addition to QCD corrections, electroweak (EW) radiative corrections 
play an important role in the W mass and width measurement; final state photon radiation 
is known to shift both quantities by O(100 MeV) [12-17]. In order to measure the W mass 
and width with high precision at a hadron collider, it is thus necessary to fully understand 
and control higher order QCD and electroweak corrections. In the last few years, significant 
progress in our understanding of the EW corrections to W boson production in hadronic 
collisions has been made. A calculation of the 0(a) EW corrections to pp — > W ± — > l^v 
(£ = e, n) in the pole approximation was presented in Ref. [18]. The complete 0(a) cor- 
rections were calculated in Ref. [19]. Two photon radiation in W production and decay 
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has been calculated in Ref. [20]. Finally, first steps towards going beyond fixed order in 
EW radiative corrections were taken in Refs. [21] and [22] by including the effects of final 
state multiphoton radiation in W production, and in Ref. [23] where final state photon ra- 
diation was added to a calculation of W boson production which includes resummed QCD 
corrections. 

In this paper we present an independent calculation of the complete O(a) EW radiative 
corrections to pp — > W ± — > i^v using the methods developed in Ref. [24], and examine 
how the non-resonant corrections neglected in Ref. [18] affect several observables of inter- 
est, in particular the measurement of the width of the W boson. Preliminary results of 
our calculation were reported in Ref. [25]. For the numerical evaluation, we use the Monte 
Carlo phase space slicing method for next-to-leading-order (NLO) calculations described in 
Ref. [26]. With the Monte Carlo method, it is easy to calculate a variety of observables si- 
multaneously and to simulate detector response. The collinear singularities associated with 
initial state photon radiation are removed by universal collinear counter terms generated 
by "renormalizing" the parton distribution functions (PDFs) [18,27-30], in complete anal- 
ogy to gluon emission in QCD. Final state charged lepton mass effects are included in our 
calculation in the following approximation. The lepton mass regularizes the collinear singu- 
larity associated with final state photon radiation. The associated mass singular logarithms 
of the form \n(s/mj), where s is the squared parton center of mass energy and mg is the 
charged lepton mass, are included in our calculation, but the very small terms of 0(mj/s) 
are neglected. 

The technical details of our calculation are described in Sec. II. The electroweak 0(a) 
corrections consist of the electroweak one-loop contributions, including virtual photons, and 
of the emission of a real photon. To regularize the ultraviolet divergences associated with 
the virtual corrections, we use dimensional regularization in the on-shell renormalization 
scheme [31]. The non-resonant corrections are small in the W pole region, but become 
important at high tv invariant masses, m(£v), due to the presence of large Sudakov-like 
electroweak logarithms of the form (a/n) In 2 (m(tv) /My) (V = W, Z) [32]. 

Numerical results for the Tevatron (pp collisions at \fs = 2 TeV) and the CERN Large 
Hadron collider (LHC, pp collisions at y/s = 14 TeV) are presented in Sec. III. For tv trans- 
verse masses above the W peak region, the non-resonant corrections reduce the differential 
cross section by 0(10%). We study in detail how these corrections affect the measurement of 
the W width from the tail of the transverse mass distribution. Using the results of Ref. [33], 
we also consider how the ratio of the W and Z cross sections, Rw/z (see Eq. (1)), and 
the transverse mass ratio of W and Z bosons are influenced. Finally, our conclusions are 
presented in Sec. IV. 

II. THE O(o) ELECTROWEAK RADIATIVE CORRECTIONS TO W 
PRODUCTION: RESONANT AND NON-RESONANT CONTRIBUTIONS 

The complete 0(a 3 ) parton level cross section of W production via the Drell-Yan mech- 
anism q^q v — > //'(t) i s gi ven by 

da^ (s, t) = da® + da virt + £ [da® F a BR + da^ 3 ] , (2) 

a = initial , f inal , 
inter f. 
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FIG. 1. The Feynman diagrams contributing to W boson production at C(a 3 ) in the Feynman 
- 't Hooft gauge (<E> + : Higgs - ghost field, u + ,u y : Faddeev-Popov-ghost fields; the non-photonic 
contribution to the W self energy insertion is symbolized by the shaded loop). An explicit repre- 
sentation of the non-photonic contribution to the W self energy insertion can be found in Ref. [34] . 

where the Born cross section, da^\ is of Breit-Wi gner form and s and t are the usual 
Mandelstam variables in the parton center of mass frame. The corresponding Feynman 
diagrams are shown in Fig. 1 and Fig. 2. da viTt describes the complete set of virtual 0(a) 
contributions consisting of the pure weak and the photonic virtual corrections: 

da virt (s, i) = ddr {0) 2Tle[F weak (s) + F 7 (s, i)} + da WZ box{s, t) . (3) 
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FIG. 2. Feynman diagrams for the W, Z box diagrams. 



Explicit expressions for the pure weak form factor F wea k and for the contribution of the 
W, Z box diagrams, d&wzbox, are given in the Appendix. The form factors F RR of Eq. (2) 
describing the initial, final state and interference contribution of real soft photon radiation, 
and the photonic form factor F 7 can be found in Ref. [24] (Appendix D). 

The soft photon region is defined by requiring that the photon energy in the parton 
center of mass frame, _E 7 , is smaller than a cutoff AE = 5 s V§/2. In this phase space region, 
the soft photon approximation can be used to calculate the cross section as long as 5 S is 
sufficiently small. Throughout the calculation the soft singularities have been regularized 
by giving the photon a fictitious mass. As usual, the unphysical photon mass dependence 
cancels in the sum of the virtual and soft photon terms, 21ZeF 1 + ^ a ^br- The IR- finite 
contribution da%^ 3 describes real photon radiation with E^ > AE. The superscript a de- 
notes the initial state, final state or interference contributions. Throughout the calculation 
of the 0(a) corrections we consider the fermions as being massless. We retain finite fermion 
masses only to regularize the collinear singularities which arise when the photon is emitted 
collinear with one of the charged fermions. Thus, d&2_, 3 and the form factors F RR and F 7 
contain large mass singular logarithms which have to be treated with special care. For final 
state photon radiation, the collinear singularity is regularized by the finite lepton mass. In 
sufficiently inclusive measurements, the mass singular logarithmic terms originating from the 
collinear singularity cancel [35]. For initial state photonic corrections, however, the mass 
singular logarithms always survive. These singularities are universal to all orders in pertur- 
bation theory and can be absorbed by a redefinition (renormalization) of the PDFs [18,27]. 
This can be done in complete analogy to the calculation of QCD radiative corrections. As 
a result, the renormalized parton distribution functions become dependent on the QED fac- 
torization scale /xqed which is controlled by the well-known Gribov-Lipatov-Altarelli-Parisi 
(GLAP) equations [36]. These universal photonic corrections can be taken into account by 
a straightforward modification [28-30] of the standard GLAP evolution equations. 

The QED induced terms in the GLAP equations lead to small corrections at the per- 
mille level to the distribution functions for most values of x and /iq ED [30,37]. Only at large 
x ~ 0.5 and large /Uq ED ~ 10 3 GeV 2 do the corrections reach the magnitude of one per cent. 

In order to treat the 0(a) initial state photonic corrections to W production in hadronic 
collisions in a consistent way, QED corrections should be incorporated in the global fitting 
of the PDFs, i.e. all data which are used to fit the parton distribution functions should be 
corrected for QED effects. Current fits [41,42] to the PDFs do not include QED corrections. 
The missing QED corrections introduce an uncertainty which, however, is probably much 
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smaller than the present experimental uncertainties on the parton distribution functions. 

Absorbing the collinear singularity into the PDFs introduces a QED factorization scheme 
dependence. The squared matrix elements for different QED factorization schemes differ by 
the finite 0(a) terms which are absorbed into the PDFs in addition to the singular terms. 
Our calculation has been carried out both in the QED MS and DIS schemes, which are 
defined analogously to the usual MS [43] and DIS [44] schemes used in QCD calculations. 
All numerical calculations in this paper are performed using the QED DIS scheme. The 
QED DIS scheme is defined by requiring the same expression for the leading and next-to- 
leading order structure function F 2 in deep inelastic scattering, which is given by the sum 
of the quark distributions. Since F 2 data are an important ingredient in extracting PDFs, 
the effect of the 0(a) QED corrections on the PDFs should be reduced in the QED DIS 
scheme. 

In the vicinity of the W resonance, the W, Z box diagrams of Fig. 2 can be neglected 
as non-resonant contributions of higher order in perturbation theory, and, as demonstrated 
in [24], a gauge invariant decomposition of the complete 0(a) contribution into a QED- 
like and a modified weak part can be performed. Unlike the Z boson case, the Feynman 
diagrams of Fig. 1 which involve a virtual photon do not represent a gauge invariant subset. 
In Ref. [24], it was demonstrated that gauge invariant contributions can be extracted from 
the infrared (IR) singular virtual photonic corrections, Fy FS (the modified Yennie-Frautschi- 
Suura form factors). These contributions can be combined with the real photon corrections 
in the soft photon region, F ER , to form gauge invariant QED-like contributions corresponding 
to initial state, final state and interference corrections. The IR finite remainder of the virtual 
photonic corrections and the pure weak one-loop corrections of Fig. 1 can be combined to 
separately gauge invariant modified weak contributions to the W boson production and 
decay processes. 

Both the QED-like and the modified weak contributions can be expressed in terms of form 
factors, Fq ED and F® eak , which multiply the Born cross section [24]. Thus, the complete 
0(a 3 ) parton level cross section of resonant W production via the Drell-Yan mechanism 
QiQi' ~^ ff'il) can be expressed in the form 



The modified weak contributions have to be evaluated at s — My^. Explicit expressions for 
the form factors Fq ED , F® eak are given in Ref. [24]. 

The non-resonant part, which was neglected in Refs. [18] and [24], can then be obtained 
from the resonant contribution, da^ 1 ^ (see Eq. (4)), and the complete 0(a) contribution, 



dfrg+U = d& (0) {1 + 2Ke(F™% 1 + F2)(M^)] 
+ ]T [da^ F^ D (s,i) + da a 2 ^]. 



(4) 



a — initial , f inal , 
inter f. 



da(°+V of Eq. (2): 



da. 



non—res 



(s,i) = da^-dai° e t 1] 



da W zbox{s, i) + da^ 2TZe 



E (Keak(s) ~ KeaM)) 



E 



FZ FS )(s,i) 



(5) 



a — initial, final, 
inter f. 
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where we have used Eq. (Al) and Fq ED = F^ R +27leF§ FS (a = initial, final, inter f.). Equa- 
tion (5) shows that the non-resonant contribution consists of the W, Z box contributions, 
the IR finite remnants of the virtual photon one-loop corrections which are not included in 
the YFS form factor, and the s-dependent parts of the pure weak vertex and self-energy 
one-loop corrections, which have been neglected when evaluating the pure weak form factor 
at s — M^r. da non _ res is free of mass singularities (terms proportional to \n(s/m 2 ), where 
nif is the mass of an initial state or final state fermion) and on-shell singularities, i.e. log- 
arithms of the form ln(]s — M^|). It represents a gauge invariant subset of the complete 
0(a) contribution to the W production process q^' — > //', and its numerical impact can 
thus be studied separately. In Fig. 3 we illustrate the effect of the non-resonant contribution 
on the total parton level cross section by showing the variation of the relative correction 
8 — & non-res / (in percent) with where is the Born cross section, 

*(o> = jl o? d - y 2 * re) 

36 1 mi) (s-M^ + (^r w y> [b) 

in the s-dependent width scheme [24]. As can be seen, the non- resonant contribution can 
be neglected in the vicinity of the W resonance but becomes increasingly important at 
large center of mass energies due to the occurrence of large, Sudakov-like logarithms of the 
form (ck/tt) In 2 (s/M^). The kink at « 200 GeV is due to the WH threshold in the W 
self-energy. 

The observable 0(a 3 ) cross section is obtained by convoluting the parton cross section 
with the quark distribution functions f q /A(x,Q 2 ) (s = xix 2 s) and summing over all quark 
flavors q and q' , 

Ms) = E/ dx lC lx 2 (f q/A ( Xl , Q 2 )f q -, /B (x 2 , Q 2 ) dM ^ + (q^ q')) (7) 

with (A,B) = (p,p) for the Tevatron and (p,p) for the LHC. The parton distribution 
functions depend on the QCD renormalization and factorization scales, \i r and /if, which we 
choose to be equal, /i r — fj,f — Q. 



III. PHENOMENOLOGICAL RESULTS 



A. Preliminaries 



We shall now discuss the phenomenological implications of the 0(a) non-resonant EW 
corrections to tv production at the Tevatron (pp collisions at y/s = 2 TeV) and the LHC 
(pp collisions at y/s = 14 TeV). For the numerical evaluation we chose the following set of 
SM input parameters [45]: 

a = 1/137.0359895, 



G u = 1.16639 x 10~ 5 GeV" 2 , 
M z = 91.1867 GeV, 
m e = 0.51099907 MeV, 
m u = 0.0464 GeV, 
m d = 0.0465 GeV, 



a s = a s (M^) = 0.121, 



m n = 0.105658389 GeV, m T = 1.777 GeV, 
m c = 1.5 GeV, m t = 174 GeV, 
m s = 0.15 GeV, m b = 4.7 GeV, 



IK 



ad | 



\V r , 



0.97 



IK. 



IK 



dc | 



0.22. 



(8) 
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FIG. 3. The relative size (in percent) of the non-resonant 0(a) corrections to the Born 
ud — > W + — ► l + v parton- level total cross section as a function of the parton center-of-mass 
energy, y/§. The parameters used are listed in Eqs. (8) - (10). 

The fermion masses only enter through loop contributions to the vector boson self energies 
and as regulators of the collinear singularities which arise in the calculation of the QED con- 
tribution. Non-zero light quark masses are only used in the calculation of the vector boson 
self energies. The light quark masses are chosen such that the value for the hadronic con- 
tribution to the photon vacuum polarization for five active flavors, Aaj^Mf) = 0.028 [4], 
which is derived from low-energy e + e~ data with the help of dispersion relations, is recovered. 
Vij are the matrix elements of the quark mixing matrix. 

The W mass and the Higgs boson mass, Mh, are related via loop corrections. A 
parametrization of the W mass which, for 10 GeV < M H < 1 TeV, deviates by at most 
0.5 MeV from the theoretical value including the full fermionic two- loop contributions is 
given in Ref . [46] . Here we use the somewhat older parametrization of Ref . [47] 

M W = K - 0.0581 in (^_) - 0.0078 In* (^_) - 0.085 (^ - l) 

_ , 18 (^_,) + ^ ((_£_)■_,) (9) 
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with = 80.3805 GeV, which was used in the analysis of the LEP data. The parametriza- 
tion of Eq. (9) reproduces the result of Ref. [47] to 0.2 MeV for 75 GeV < M H < 350 GeV. 
For the numerical discussion we choose 



M H = 120 GeV, 



(10) 



which is consistent with current direct [2] and indirect bounds [1,5], and work in the s- 
dependent width scheme. For the input parameters listed in Eq. (8) we obtain M\y = 
80.3612 GeV. Higher order (irreducible) corrections connected with the p-parameter are 
taken into account in our calculation by replacing 



8M% 
"mJ" 



M 2 W 



5M 2 Z 



&Mw A n HO 



(11) 



in the renormalized Z and 7, Z self-energies as described in Appendix B of Ref. [18]. The W 
width and the — ■> ff'i'j) amplitude are calculated in the G u scheme. In the scheme, 
the fine structure constant, a, is replaced with 



a 



V2G^M } 



2 

w 



1 - 



w 



7Y 



(12) 



and M w is obtained from Eq. (9). For the W width, including 0(a) EW (see also Ref. [38]) 
and QCD corrections up to 0(a\) [18,39,40], we obtain Y w = 2.0721 GeV. 

In the numerical results presented below, we also take into account the leading 0(a 2 ) 
weak corrections to resonant W production which can be obtained by performing the fol- 
lowing replacement in Eq. (4): 



2TZe{Fl 



weak 



p final} 1 ;i f'2 



weak 



){M* W ) 



1 ^ 
2 



I- f? final / n*-2 
~ r weak \ m \ 



W . 



(13) 



The dominant non-photonic electroweak corrections to resonant W production can be 
taken into account in the effective Born approximation (EBA) where the W cross section, 
a EBA , is calculated by replacing a in Eq. (6) with the expression given in Eq. (12) and using 
the W mass obtained from Eq. (9). In the following, we shall use the EBA cross section as 
a reference when discussing how the 0(a) EW corrections affect physical observables. 

Following a brief comparison with the calculation of Ref. [19], we discuss the impact 
of the non-resonant EW corrections on the transverse mass distribution and how they af- 
fect the value of the H^-width extracted from that distribution. We then consider the W 
production cross section, the W to Z cross section ratio and the W to Z transverse mass 
ratio. To compute the hadronic cross section we use the MRSR2 set of parton distribution 
functions [48], and take the renormalization scale, fj, r , and the QED and QCD factoriza- 
tion scales, /xqed and //qcd, to be [/% = /!q ED = Aiq CD = M^. The detector acceptance 
is simulated by imposing the following transverse momentum (p T ) and pseudo-rapidity (77) 
cuts: 



p T (£) > 20 GeV, 



1 77 



< 2.5, 



i = e, fi, 



(14) 



j> T > 20 GeV, 



(15) 
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where f T is the missing transverse momentum originating from the neutrino. These cuts 
approximately model the acceptance of the CDF II [49] and D0 [50] detectors at the Teva- 
tron, and the ATLAS [51] and CMS [52] detectors at the LHC. Uncertainties in the energy 
measurements of the charged leptons in the detector are simulated in the calculation by 
Gaussian smearing of the particle four-momentum vector with standard deviation a which 
depends on the particle type and the detector. The numerical results presented here were 
calculated using a values based on the D0 and ATLAS specifications. 

The granularity of the detectors and the size of the electromagnetic showers in the 
calorimeter make it difficult to discriminate between electrons and photons with a small 
opening angle. In such cases we recombine the four-momentum vectors of the electron and 
photon to an effective electron four-momentum vector. The exact recombination procedure 
is detector dependent. For calculations performed at Tevatron energies we use a procedure 
similar to that used by the D0 Collaboration in Run I, requiring that the electron and 
photon momentum four-vectors are combined into an effective electron momentum four- 
vector if their separation in the pseudorapidity - azimuthal angle plane, 

AR(e, 7 ) = v / (A7 ? ( e ,7)) 2 + (A0(e ) 7)) 2 , (16) 

is Ai2(e,7) < 0.2. For 0.2 < AR(e,j) < 0.4 events are rejected if E y > 0.15 E e . Here 
E 1 (E e ) is the energy of the photon (electron) in the laboratory frame. For events with 
0.2 < AR(e, 7) < 0.3 and E 1 < 0.15 E e , the electron and photon momentum four- vectors are 
again combined. At LHC energies, we recombine the electron and photon four-momentum 
vectors if Ai?(e, 7) < 0.07, similar to the resolution expected for ATLAS [51]. Recombining 
the electron and photon four-momentum vectors eliminates the mass singular logarithmic 
terms originating from final state photon radiation and strongly reduces the size of the 
QED-like final state corrections [18]. 

Muons are identified by hits in the muon chambers and the requirement that the associ- 
ated track is consistent with a minimum ionizing particle. This limits the photon energy for 
small muon - photon opening angles. For muons at the Tevatron, we again adopt the D0 
specifications and require that the energy of the photon is E 1 < 2 GeV for Ai2(/i, 7) < 0.2, 
and E y < 6 GeV for 0.2 < Aify/,7) < 0.6. At the LHC, following Ref. [51], we require the 
photon energy to be smaller than EJ = 5 GeV if AR(/j,, 7) < 0.3. The cut on the photon 
energy increases the size of the photonic corrections for m(/j,u) > 100 GeV [18]. For future 
reference, we summarize the lepton identification requirements in Table I. 

We use the input parameters listed in Eq. (8) and impose the cuts and lepton identifica- 
tion requirements described above in all subsequent numerical simulations, unless explicitly 
noted otherwise. 

B. Comparison with Ref. [19] 

As mentioned before, the matrix elements of the full 0(a) EW corrections to qq' — > lv+X 
were presented in Ref. [19], together with a discussion of how EW radiative corrections 
influence the pr distribution of the charged lepton, and the transverse mass distribution. In 
this section we compare the integrated cross sections given in Ref. [19] for different ranges 
in pr(£) at the Tevatron and LHC with the results of our calculation, using the input 
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TABLE I. Summary of lepton identification requirements. 



Tevatron 




electrons 


muons 


combine e and 7 momentum four vectors if 


reject events with E 7 > 2 GeV 


A J R(e,7) < 0.2 and if E y < 0.15 E e for 0.2 < AR(e,~/) < 0.3 


for Ai?(/x,7) < 0.2 


reject events with Ey > 0.15 E e 


reject events with > 6 GeV 


for 0.2 < AR(e,j) < 0.4 


for 0.2 < &R(ji,i) < 0.6 


LHC 




electrons 


muons 


combine e and 7 momentum four vectors if 


reject events with E y > 5 GeV 


Ai?(e, 7 ) < 0.07 


for AJ?(/x,7) < 0.3 



parameters, PDFs, cuts, and the lepton - photon recombination procedure of Ref. [19]. The 
results are shown in Table II. The table lists the lowest order cross section in the scheme, 
and the relative corrections, 



a O(a 3 ) _ (0) 

— ^T- < 17 > 



for the full O(o) EW corrections, and 



a PA - <j(°) 

for the 0(a) EW corrections in the pole approximation. Our results are found to agree within 
the statistical accuracy of the Monte Carlo integration with those obtained in Ref. [19] over 
the entire lepton p T range. 



C. Non-resonant corrections to the transverse mass distribution 

Table II shows that the non-resonant EW corrections quickly become important for 
transverse momenta above the Jacobian peak region, pr(£) > M w /2 ps 40 GeV. The devia- 
tion of the 0(a 3 ) cross section in the pole approximation from the full NLO result at large 
Pt(£) is due to large Sudakov-like electroweak logarithms of the form (en /it) In (m(iv)jMy) 
(V = W, Z) [32] which mostly arise from the contribution of the W, Z box diagrams and to 
a lesser extend from the energy dependence of the form factors F£ eak and F 7 — J2 a ^YFS- 
The same qualitative behavior is also expected for the transverse mass distribution, which 
is used to determine the W mass and width [14]. The transverse mass is defined by 

M T = ^2 Pt (£) Pt (u)(1-cos^) , (19) 

where pr(v) is the transverse momentum of the neutrino, and (f> eu is the angle between the 
charged lepton and the neutrino in the transverse plane. The neutrino transverse momentum 
is identified with the missing transverse momentum, $ T , in the event. 
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TABLE II. Integrated lowest order cross sections in the scheme, the relative corrections 
for the pole approximation, Spa, and the full 0{a) EW corrections, S, for several ranges of Pt(£)- 
Shown are the results of our calculation and of Ref. [19] for a) pp — ► iv at \fs = 2 TeV and b) 
pp — > £^ at \/s = 14 TeV. The statistical uncertainties of the Monte Carlo integration are also 
shown. The lepton and photon momenta are recombined for small values of AR(£, 7) using the 
simplified procedure described in Ref. [19]. We also use the same input parameters, PDFs and cuts 
as Ref. [19]. 



a) pp -> Iv, y/s = 2 TeV 



p T {£) (GeV) 


25-oo 


50-oo 


75 — 00 


100 - 00 


200 - 00 


300- 00 


crW (pb) Ref. [19] 


407.03(5) 


2.481(1) 


0.3991(1) 


0.1305(1) 


0.006020(2) 


0.0004821(1) 


o-(°) (pb) this calc. 


407.02(7) 


2.4817(6) 


0.39926(9) 


0.13058(3) 


0.006017(2) 


0.0004821(3) 


5 (%) Ref. [19] 


-1.8(1) 


-2.7(1) 


-4.8(1) 


-6.3(1) 


-10.4(1) 


-13.6(1) 


S (%) this calc. 


-1.7(1) 


-2.5(1) 


-4.7(1) 


-6.1(1) 


-10.1(1) 


-13.3(1) 


5p A (%) Ref. [19] 


-1.7(1) 


-1.6(1) 


-2.3(1) 


-2.5(1) 


-3.3(1) 


-3.9(1) 


$pa (%) this calc. 


-1.7(1) 


-1.5(1) 


-2.2(1) 


-2.4(1) 


-3.1(1) 


-3.7(1) 






b) pp - 


+ £u : ^ = 14 TeV 






pr(^) (GeV) 


25-oo 


50-00 


100 - 00 


200 - 00 


500 - 00 


1000 - 00 


(pb) Ref. [19] 


1933.5(3) 


11.50(1) 


0.8198(4) 


0.1015(1) 


0.005277(1) 


0.0003019(1) 


cj( ) (pb) this calc. 


1933.4(3) 


11.499(2) 


0.8202(1) 


0.10155(2) 


0.005277(1) 


0.0003019(1) 


<5 (%) Ref. [19] 


-1.8(1) 


-2.7(1) 


-6.2(1) 


-10.2(1) 


-19.6(1) 


-29.6(1) 


<5 (%) this calc. 


-1.8(1) 


-2.3(1) 


-6.0(1) 


-10.1(1) 


-19.1(1) 


-28.6(1) 


5p A (%) Ref. [19] 


-1.8(1) 


-1.5(1) 


-1.6(1) 


-1.6(1) 


-2.4(1) 


-3.6(1) 


<W (%) this calc. 


-1.8(1) 


-1.2(1) 


-1.5(1) 


-1.6(1) 


-2.4(1) 


-3.4(1) 
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FIG. 4. The ratio [da°^ /dM T ]/[da EBA /dM T ] as a function of the transverse mass for a) 
pp — ► e + v e ('f) and b) pp — > fi + iy^(^f) at = 2 TeV. The solid lines show the ratio of the complete 
0(o?) electroweak and the EBA differential cross section. The dashed lines display the correspond- 
ing ratio for the case where only the resonant 0(a) EW corrections (see Eq. (4)) are taken into 
account (pole approximation). The dotted lines show the ratio when the pp — > W + (-^ £v)Z(—> vv) 
background is included in addition to the complete O(a) EW corrections. The cuts and lepton 
identification requirements imposed are described in Sec. Ill A. 

The ratio of the complete 0(a 3 ) electroweak and the EBA differential cross section as 
a function of is shown in Fig. 4. In order to make the effect of the non-resonant weak 
corrections more transparent, we also show the corresponding ratio for the case of the 0(a) 
EW corrections in the pole approximation (dashed lines) [18]. For < Mw, the pole 
approximation is seen to very well represent the complete electroweak 0(a) corrections. In 
this region, the shape change in the Mt distribution is largely due to the contribution of 
the final state QED-like corrections. Due to the recombination of electrons and photons, 
the EW corrections in the pole approximation reduce the ev(pj) differential cross section by 
only 2 — 3% over the transverse mass region considered. In the muon case, the cut on the 
photon energy for photons which have a small opening angle with the muon reduces the 
hard photon part of the 0(a 3 ) ^(7) cross section. As a result, the QED-like corrections 
are much more pronounced and display a much stronger dependence on the transverse mass 
than in the electron case. Without taking the lepton identification criteria of Table I into 
account, the QED-like corrections in the electron case are larger, due to the mass singular 
logarithmic corrections which originate from final state photon radiation. The non-resonant 
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FIG. 5. The ratio [da°( aS ^ / ' dMT\/[da EBA / 'cIMt] as a function of the transverse mass for a) 
pp — > e + v e ( / y) and b) pp — ► p + v^("f) at -^/s = 14 TeV. The solid lines show the ratio of the complete 
C(a 3 ) electroweak and the EBA differential cross section. The dashed lines display the correspond- 
ing ratio for the case where only the resonant 0(a) EW corrections (see Eq. (4)) are taken into ac- 
count (pole approximation). The dotted lines show the ratio when the pp — ► W + {-^ i + u)Z{-^ vv) 
background is included in addition to the complete 0(a 3 ) EW corrections. The cuts and lepton 
identification requirements imposed are described in Sec. Ill A. 

EW corrections are seen to increase rapidly in size with M T . For M T = 300 GeV, they 
reduce the cross section by about 4% at the Tevatron. The dotted lines, finally, show cross 
section ratio when the pp — > W + {~^ i + v)Z(-^ vv) background is included in addition to 
the complete 0(a) EW corrections. The W Z background is seen to be much smaller than 
the 0(a) EW corrections. 

Figure 5 shows the corresponding results for the LHC in the high transverse mass tail. 
For Mt > 1.5 TeV, the non-resonant 0(a) EW corrections reduce the differential cross 
section by 20% or more and thus are of the same size as the 0(a s ) corrections. This is 
larger than the expected statistical uncertainty in a 200 GeV bin centered at Mt = 1.5 TeV 
for 100 fb _1 . It will thus be important to take into account the non-resonant weak correc- 
tions when searching for new heavy W bosons, such as Kaluza-Klein excitations of the W 
appearing in TeV-scale models with extra dimensions [53], at the LHC. The results shown 
in Fig. 5, however, should be interpreted with caution. Since the non-resonant weak correc- 
tions become large for transverse masses above 1 TeV, they need to be resummed in order 
to obtain accurate predictions in this phase space region (for a recent review of the resum- 
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mation of electroweak Sudakov-like logarithms see Ref. [54]). Although the resummation of 
electroweak Sudakov-like logarithms in general four fermion electroweak processes has been 
discussed in the literature [55], a calculation of Iv production in hadronic collisions which 
includes resummation of electroweak logarithms has not been carried out yet. 



D. Non-resonant EW radiative corrections and the W width 

We now discuss how the non-resonant 0(a) EW corrections affect the W width deter- 
mined from the high transverse mass tail. In Run I of the Tevatron, Yyy has been measured 
by the CDF and D0 Collaborations using this technique with a combined uncertainty of 
105 MeV [14]. In Run II, with an integrated luminosity of 2 fb" 1 , one expects to achieve a 
precision of 50 MeV per lepton channel and experiment [7]. Assuming that the error corre- 
lation of CDF and D0 data, as in Run I, is small [14] for the W^-width measurement, this 
results in an expected overall uncertainty of STw = 25 — 30 MeV. In the experimental anal- 
ysis, the measured transverse mass distribution is compared with the theoretical prediction 
for various values of the W width where the total cross section has been normalized to the 
experimental value [16,17]. This is equivalent to analyzing the normalized M T distribution, 

where Otot is the total iv(pj) cross section within cuts. At lowest order, the W width enters 
the cross section in the form of the squared W propagator (see Eq. (6)), 

\D w (s)\ 2 = [(§ - M 2 W ) 2 + s 2 T 2 w /M 2 w }-\ (21) 

The lowest order total tv production cross section thus is proportional to \ As a result, 
da/dM T is proportional to T w for M T 3> M w . This is clearly displayed in the ratio of the 
normalized M T distribution for arbitrary T w to the normalized M T distribution in the SM, 



n{v w ) = 




where = 2.072 GeV is the SM W width (see Sec. Ill A). K{T W ) in the pole 

approximation, TZpa(^w), is shown in Fig. 6 for pp — > e + u('~f) at the Tevatron with 
IV = Tffi - 10 MeV = 2.062 GeV (dotted line) and T w = T^ 1 - 30 MeV = 2.042 GeV 
(dashed line). Qualitatively similar results are obtained for the muon final state and at LHC 
energies. 1Z PA (T W ) is seen to be almost constant and approximately equal to T w /T^ for 
Mt > 100 GeV. For < Mw, the ratio is almost independent of T^, and very close 
to one. If no detector resolution effects are taken into account, dajdM^ is proportional to 
l/TV for Mt = M\y at lowest order, resulting in an enhancement proportional to T^ 1 /Tw 
in lZ{Tw)- However, detector resolution effects largely dilute this effect. Above the W 
mass, 1Z P a(Xw) rapidly drops from 1Z P a(Xw) ~ 1 to TZpa(^w) ~ ^w/^w 1 ■ The M T range 
over which the transition occurs sensitively depends on the missing transverse momentum 
resolution of the detector. 
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FIG. 6. The ratio TZiTw) (see text for definition) as a function of My for PP — ► e + v(^) at 
the Tevatron. Results are shown for T w = T^ 1 - 10 MeV = 2.062 GeV (dotted line) and 
IV = ^w 1 ~ 30 MeV = 2.042 GeV (dashed line) with the cross sections calculated in the pole 



approximation. The solid line displays the ratio TZ P ^ = [da°^ a ) / 'dMT\/[da PA / 'cIMt] in the SM. 
The data points and error bars indicate the measurements and statistical uncertainties expected 
for 20 GeV bins and an integrated luminosity of 2 fb _1 , assuming that the data are 



for n 



0{a A 
PA 



described by the SM prediction including the full 0(a) EW corrections. The cuts and lepton 
identification requirements imposed are described in Sec. Ill A. 



The solid line in Fig. 6 shows the ratio 



0{a 3 ) 
PA 



dM T 
da PA 



(23) 



for Y w = Ty^ 1 , where da 0{pL ) /dM T and da PA / dM T are the SM normalized M T distributions 
for the complete 0(a 3 ) calculation and in the pole approximation, respectively. The data 
points and error bars in Fig. 6 indicate the measurements and statistical uncertainties ex- 
pected for IZ^a ^ f° r 20 GeV bins and an integrated luminosity of 2 fb" 1 , assuming that the 
data are described by the SM prediction including the full 0(a) EW corrections. Due to the 
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non-resonant EW corrections, TZpX gradually decreases for M T > 90 GeV, the region used 
by the Tevatron experiments to extract the W width from Run I data [16,17]. Although 
the shape change of the transverse mass distribution due to the non-resonant 0(a) EW 
correction differs significantly from that caused by a non-standard W width, the expected 
statistical uncertainties make it difficult to distinguish between a small negative shift in IV 
and the effect of non-resonant 0(a) EW corrections. A \ 2 analysis shows that, were the 
non- resonant 0(a) EW corrections ignored, the value of the W width extracted from the 
M T > 90 GeV data region would be shifted by 

AIV = -7.2 MeV. (24) 

Since the M T distribution depends little on the detector resolution for M T > 90 GeV, the 
shift in T w is almost independent of these effects [7]. For the precision of T w expected in 
Run II, a difference of ^ 7 MeV in the extracted value of the W width cannot be ignored, 
and the complete 0(a 3 ) calculation should be used to compare theory and data. 



E. Non-resonant corrections to the W boson cross section and the W to Z cross 

section ratio 

As mentioned in the Introduction, the W width can also be determined from the cross 
section ratio Rw/z (see Eq. (1)), together with the theoretical prediction for the ratio of the 
total W and Z production cross sections, the LEP measurement of the Z — > £ + £~ branching 
ratio, and the SM prediction for the W — > tv decay width. Since the QCD corrections to 
W and Z production are very similar, they cancel almost perfectly in the W to Z cross 
section ratio; the 0(a s ) corrections to Rw/z are of 0(1%) or less, depending on the set of 
parton distribution functions used [56]. In addition many experimental uncertainties, such 
as the luminosity uncertainty, cancel in the cross section ratio. Accurate knowledge of how 
electroweak corrections affect Rw/z is thus very important. 

The W cross section may be used as a luminosity monitor in the future [57,58]. This 
requires that the W cross section is reliably computed with small uncertainty. In order to 
achieve this, it is essential to know how EW radiative corrections affect the W — > Iv cross 
section. 

The size of the 0(a) electroweak corrections to the total pp — > ivX cross section and to 
Rw/z is sensitive to the acceptance cuts and whether lepton identification requirements are 
taken into account or not. In Table III, we list the electroweak i^-factor, 

a EBA( pP _ w tv) ' y °) 

and the correction factor for Rw/z, 

r>0(a 3 ) 

R EW _ U W/Z (26) 
^W/Z 

for the acceptance cuts listed in Eqs. (14) and (15) with and without taking the lepton 
identification requirements of Table I into account. For K EW , we also list the result obtained 
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TABLE III. The electroweak K-factor K EW = a ^ (pp -► W -► £vX)/a EBA {pp W ^ £v) 
(£ = e, n) and the correction factor to i^/z, = R w"z l R wjz > with 75 GeV < m (^~) 

< 105 GeV, for pp collisions at y/s = 2 TeV. Shown are the predictions without and with the lepton 
identification requirements of Table I taken into account. The cuts imposed are listed in Eqs. (14) 
and (15). The energy and momentum resolutions used are described in Sec. Ill A. 





without lepton id. 


with lepton id. 




requirements 


requirements 


k ew ( pp _^ e + uX ), full 0(a 3 ) 


0.963 


0.984 


K EW {pp — > e + uX), pole appr. 


0.963 


0.984 


k ew ( pp _^ ^uX), full O(o 3 ) 


0.979 


0.944 


K EW {pp — > /j, + i/X), pole appr. 


0.978 


0.943 


K% w (e), full 0(a 3 ) 


1.024 


0.992 


K^ w (ji), full 0(a 3 ) 


1.002 


1.045 



in the pole approximation. To compute the 0{a 3 ) Z boson cross section entering Rw/z, 
we use the full 0{a 3 ) calculation of di-lepton production in hadronic collisions described in 
Ref. [33]. We include photon exchange and 7Z interference effects, and impose a cut on the 
di-lepton invariant mass of 75 GeV < m{£ + £~) < 105 GeV. The definition of the effective 
Born approximation in the Z case is given in Ref. [33]. The values listed in Table III update 
the results presented in Ref. [18] which were obtained using the pole approximation for the 
W — > iv cross section and only included QED corrections in the Z boson case. 

From the results listed in Table III one observes that the non-resonant O{o) EW cor- 
rections have a very small effect on the W boson cross section; they change the electroweak 
K factor by (9(1CT 3 ) or less. The full 0(a) EW corrections decrease the W cross section 
and increase Rw/z by several per cent for the cuts imposed. When lepton identification 
requirements are included, the corrections are reduced in the electron case and enhanced 
in the muon case. Unlike the QCD corrections, the electroweak corrections do not cancel 
in Rw/z- In Z —> £ + £~ both leptons can emit photons, whereas only the charged lepton 
radiates in W — > Iv decays. Since final state photonic corrections are the dominating con- 
tribution to the 0{a) EW corrections, the 0{a) corrections to the W and Z cross sections 
are quite different, and therefore do not cancel in Rw/z- They are of the same size as the 
QCD corrections to Rw/z- The size of the 0(a) EW corrections to the total pp — > £vX 
cross section and to Rw/z is similar to the statistical uncertainty from 72 pb -1 of data from 
Run II [59]. It is thus important to take the EW radiative corrections into account in the 
Run II data analysis. 

F. The full 0{a) electroweak corrections to the W to Z transverse mass ratio 

Since detectors cannot directly detect the neutrinos produced in the leptonic W boson 
decays, W — > £is, and cannot measure the longitudinal component of the recoil momentum, 
there is insufficient information to reconstruct the invariant mass of the W boson. Instead, 
the transverse mass distribution of the final state lepton pair, or the transverse momentum 
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distribution of the charged lepton are used [12,13] to extract Mw- The Mt distribution has 
the advantage of being invariant under transverse boosts to first order in the velocity of the 
W boson. On the other hand, the transverse mass depends on an accurate reconstruction of 
the neutrino direction which leads to a set of experimental requirements which are difficult 
in practice to control [7]. 

The W mass can also be determined from the ratio of the transverse mass distributions of 
the W and Z boson [60-62]. The advantage of this method is that one can cancel common 
scale factors in ratios and directly determine M w /M z , which can be compared with the 
precise value of M z from the LEP experiments. The downside of the ratio method is that 
the statistical precision of the Z sample is directly propagated into the resultant overall 
uncertainty of Mw- At high luminosities, or when a detailed understanding of the detector 
response is not available, the transverse mass ratio of W to Z bosons offers advantages in 
determining the W mass over the M T distribution [7]. 

The transverse mass ratio of W and Z bosons is defined as 



Aw(X M — X 



Mr, 



Rm t (X Mt ) - , (v z T _ v — V ' ( 27 ) 



where A v (V = W, Z) is the differential cross section 



MX V M T ) = ^r- (28) 



with respect to the scaled transverse mass, 

Xl T - % . (29, 

The transverse mass of the lepton pair in Z boson events is defined in complete analogy to 
Eq. (19): 

M# = v /2p T (£+)p T (£-)(l-cos0«) , (30) 

where 4> u is the angle between the two charged leptons in the transverse plane. The 0(a) EW 
radiative corrections to Rm t were calculated in Ref. [18] in the approximation where only 
the QED corrections were taken into account for Z — > £ + £~ and the pole approximation was 
used for W — > Iv. Here we present results which include the complete 0(a) EW radiative 
corrections. 

The ratio of the 0(a 3 ) and the W to Z transverse mass ratio in the effective Born 
approximation is shown in Fig. 7. To calculate the 0(a) electroweak corrections to Z 
boson production, we again use the results of Ref. [33]. As before, photon exchange and 
7Z interference effects are included and an additional cut on the di-lepton invariant mass 
of 75 GeV < m(£ + £~) < 105 GeV has been imposed. The complete 0(a) EW radiative 
corrections (solid lines) are uniformly about 1% smaller than those obtained when the pole 
approximation is used for W production and only QED corrections to Z boson production 
are taken into account (dashed lines). Most of this effect originates from the genuine weak 
corrections to the Z boson cross section. In the electron case, the 0(a) EW corrections 
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FIG. 7. Ratio of the C(a 3 ) and the W + to Z transverse mass ratio in the effective Born 
approximation as a function of the scaled transverse mass, Xm t , for pp collisions at yfs = 2 TeV, 
a) for electron, and b) for muon final states. The solid lines show the result when the complete 
O(a) EW radiative corrections are taken into account. The dashed lines give the result for the 
EW radiative corrections to W production in the pole approximation, and when only the QED 
corrections to Z boson production are taken into account (as in Ref. [18]). The cuts and lepton 
identification requirements imposed are described in Sec. Ill A. For pp — ► £ + l~(^f), we in addition 
require the di- lepton invariant mass to satisfy the constraint 75 GeV < m(l + l~) < 105 GeV. 

change Rm t by 1 — 2%. They are larger in the muon case (3 — 7%) due to the lepton 
identification requirements which reduce the hard photon part of the 0(a 3 ) cross section 
and thus enhance the effect of the virtual corrections. In the resonance region, Xm t ~ 1, 
hard photonic corrections reduce both the W and Z boson cross sections. As mentioned 
before, in pp — > iv{^{) only one of the two leptons can emit a photon, whereas both leptons 
in pp — > £ + £~('~f) can radiate. As a result, the effect of the hard photonic corrections is more 
pronounced in the Z case, resulting in a resonance- like enhancement of Rm t for Xm t ~ 1- 



IV. CONCLUSIONS 

We have presented a calculation of the O(a) corrections to pp — » — > l~^v based 
on the complete set of one-loop Feynman diagrams contributing to Iv production, using 
the methods developed in Ref. [24]. The calculation is based on a combination of analytic 
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and Monte Carlo integration techniques. Lepton mass effects are included in the approxi- 
mation where only mass singular terms originating from the collinear singularity associated 
with final state photon radiation are retained. The ultraviolet divergences associated with 
the virtual corrections are regularized using dimensional regularization and the on-shell 
renormalization scheme [31]. The cross sections obtained using our matrix elements were 
found to be in very good agreement with those obtained in Ref. [19]. 

Since the structure of the full 0(a 3 ) matrix elements and the 0(a) EW corrections in 
the pole approximation were discussed in detail in earlier papers [18,19], we concentrated on 
the phenomenological effects of the non-resonant corrections. The non-resonant corrections 
were found to have a very small effect on the total W cross section and the transverse 
mass distribution in the region Mt < My/- However, they increase rapidly in magnitude 
with Mt above the W peak, due to the presence of Sudakov-like electroweak logarithms. 
Although these corrections are of moderate size for transverse masses accessible at the 
Tevatron, they induce a shift of AIV ~ — 7 MeV in the W width extracted from the tail of 
the transverse mass distribution. Comparison with the expected overall precision of about 
STw = 25 — 30 MeV in Run II of the Tevatron shows that it will be necessary to take 
the non-resonant electroweak corrections into account in the data analysis. For transverse 
masses in the TeV region which play an important role in new physics searches at the LHC, 
the non- resonant electroweak corrections are of the same size as the 0(a s ) corrections. The 
strong increase of these corrections with M T requires that they are resummed. No such 
calculation exists yet for iv production in hadronic collisions. 

We also updated the results of the electroweak i^-factor for the W/Z cross section ratio 
Rw/z an d the W to Z transverse mass ratio given in Ref. [18] to include the complete 0(a) 
EW radiative corrections and the 0(g A ml j 'My/) corrections to the Z boson cross section [33]. 
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The pure weak vertex and self energy one-loop corrections to the process described by 
the form factor F weak of Eq. (2) can be decomposed into initial and final state contributions 
as follows 
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APPENDIX A: THE PURELY WEAK CONTRIBUTION 




rinitial 



i final 



(Al) 



with the final state contribution 
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1 spW,weak,~\ _ spW,weak , ** 2 \ i 
j=I,II,III A S - M w 



(A2) 



Performing the substitution (/, /') — > (i, i') yields the corresponding initial state form factor 
F™ea k al (s). The explicit expressions for the various terms in Eq. (A2) can be found in 
Ref. [24]. 

The pure weak box contribution consisting of the WZ box diagrams of Fig. 2 in the limit 
of massless external fermions reads 



da WZbox (s, t) = - dM°\s, t)(s - M 2 W ) Ke(5B vlV2 + SB^f) (A3) 



V1,V2=Z,W 

with 



5B vlV2 (s, t) = \ V1V2 [2D° 2 + t(D\ + D\ + D\ + D\ + Df + D\ 2 ) 
5B™ s 2 ed (s, u) = \ C V1V2 \8D° 2 + u(Dl + 2D\ + D\ + 2{D\ + D? + D 12 )) - 2sD 



13 > 



(A4) 



where L>/ (fc) = £>f (fe) (s, t, 0, M V i, 0, M V2 ) in case of the box, 5B V1V2 , and L>/ (fc) = 

in case of the crossed box contribution, 5By^y 2 d . \v1v2 and 
X V iV2 contain the dependence on the vector and axial vector parts of the fermion-Z boson 
couplings which, in the notation of Fig. 8, are given by 

\ wz = 2(vi + ai)(v 3 + a 3 ), 
Xzw = 2(w 2 + a 2 )(v 4 + a 4 ), 
\ c wz = 2(vi + ai)(u 4 + a 4 ), 

A^ = 2(v 2 + a 2 )(t»3 + a 3 ), (A5) 

with Vi = (I\ — 2s 2 Qi) I '{2s w c w ) and Oj = /3/(2s w c w ) parametrizing the Z// (/ = £, g) 
couplings. Here, Qi and /g denote the charge and third component of the weak isospin 
quantum numbers of the fermion i, and s w = sin9 w , c w = cos9 w with 9 W being the weak 
mixing angle. 

The box integrals arising in the calculation of the Feynman diagrams of Fig. 8 are of the 
form 



d A k (lMMk v ) 



16tt 2V u ' ' ' J (2tt) 4 [P + ie][(k - p 2 ) 2 - M vi }[{k + p k f + ie][(k + Pl ) 2 - M V2 



2 

(AC) 



The explicit decomposition of the vectorial and tensorial four point functions 

D" = -£D\ + vlD\ + p?£>? 

= f 2 f 2 D\ + V lf k D\ + pf^Df + gTDl 

- (p&* + lfa®D? - (rfrt + rt&)D? + {fkVl + P1PI)D? (A7) 
can be found in Ref. [63]. 
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FIG. 8. Notations used in the calculation of the VI, V2 box diagrams. 
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